Exact crossover Green function in the two-channel and two-impurity Kondo models 
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Symmetry-breaking perturbations destabilize the critical points of the two-channel and two- 
impurity Kondo models, thereby leading to a crossover from non-Fermi liquid behavior to stan- 
dard Fermi liquid physics. Here we use an analogy between this crossover and one occurring in 
the boundary Ising model to calculate the full crossover Green function analytically. In remarkable 
agreement with our numerical renormalization group calculations, the single exact function applies 
for an arbitrary mixture of the relevant perturbations in each model. This rich behavior resulting 
from finite channel asymmetry, inter- lead charge transfer and/or magnetic field should be observable 
in quantum dot or tunneling experiments. 
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The most basic quantum impurity model exhibiting 
non-Fermi liquid (NFL) behavior is arguably the two- 
channel Kondo (2CK) model [l| , describing the symmet- 
ric antiferromagnetic coupling of a local spin-^ impurity 
to two equivalent but independent conduction channels. 
The resulting ground state possesses various intriguing 
properties, including notably a residual entropy [2j of 
ifeg ln(2) and conductance that approaches its T = 
value as y/T (for a review see Ref. 0). 

The same behaviour is predicted at the critical point of 
the two-impurity Kondo (2IK) model Q|. The tendency 
to form a trivial local singlet state is favoured by an ex- 
change coupling acting directly between the impurities; 
while the coupling of each impurity to its own metal- 
lic lead favours separate single-channel Kondo screening. 
The resulting competition gives rise to a critical point 
that is closely related to the 2CK state. 

The central difficulty in realizing experimentally the 
NFL physics of either model is the extreme delicacy of 
the 2CK fixed point (FP) to various symmetry-breaking 
perturbations. Channel asymmetry, magnetic field and 
inter-lead charge transfer processes all destabilise the 
2CK FP and destroy NFL behaviour in both 2CK and 
2IK models. 

Tremendous efforts have been made to suppress these 
relevant perturbations in order to observe the character- 
istic NFL behaviour in a real 2CK device. The quan- 
tum dot system realized recently in Ref. [5j shows unam- 
biguous signatures associated with flow to the 2CK FP. 
But in any real system, the presence of destabilising per- 
turbations is totally inevitable, leading ultimately to a 
crossover from NFL behaviour to standard Fermi liquid 
physics on the lowest energy scales. 

In this Letter we demonstrate that an arbitrary mix- 
ture of the relevant perturbations in either the 2CK or 
2IK model leads to low-energy behaviour of the impu- 
rity Green function that is wholly characteristic of the 
incipient 2CK state, and which can not be extracted 
from the Bethe-ansatz solution (i|. We derive a single 
exact function to describe this crossover, which agrees 
perfectly with our full numerical renormalization group 



(NRG) calculations for both models, and whose rich be- 
haviour should be directly observable in quantum dot or 
tunneling experiments. 

Models and perturbations.- We consider the standard 
2CK and 2IK models, 

H 2 ck = H Q +JS ■ (s L + s r) + 8H 2 ck, (1) 
H 2 ik = H Q +J(S L ■ s L + S R ■ s 0R ) + KS L ■ S R + SH 2IK , 

(2) 



where H, 



ii 



describes two free con- 
duction electron channels a = L/R, with spin density 

SOa = E CT(7 ' ^l™ (k^'Wocr'a ^0™ = Efe 

coupled to one spin-^ impurity S (2CK) or two impurity 

spins Sl,r (2IK). For 8H 2 ck — 0, the ground state of 
H 2 ck is described by the 2CK FP. Likewise, a critical 
inter-impurity coupling K c can be found such that the 
ground state of R 2 ik is similarly described by the 2CK 
FP for SH 2IK = 0. 

Relevant perturbations to each model have been iden- 
tified from conformal field theory (CFT) 0, 0| and are 
generically present in experiment. Specifically, 
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(3) 

describes L/R channel asymmetry in the 2CK model for 
A z 0, while charge transfer between the leads is em- 
bodied in the A x and A y components of the first term 
[here f (a) are the Pauli matrices in the channel (spin) 
sector]. The second term describes a magnetic field act- 
ing on the impurity. For the 2IK model, the critical point 
is destabilized by finite (K c — K), and also through 
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{V lr ^ l ^r+H..c.)+B s -{S l -Sr), (4) 



where Vlr describes electron tunneling between the leads 
and B s the application of a staggered magnetic field. 

These perturbations generate a new energy scale, T*, 
characterizing the flow away from the 2CK FP, and 
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Figure 1: Spectrum i CT (6j) vs. uj /Tk for the 2CK model (upper panels) and the 2IK model (lower panels) at T = in the presence 
of various perturbations. Entire frequency dependence calculated by NRG (black lines) ; low-energy uj <C Tk behaviour in each 
case compared with exact crossover function Eq. (red lines). All results presented for v. J — 0.25. Left: effect of channel 
asymmetry A z / (2CK) or deviation from critical coupling (K c - K ) (2IK). Specifically, AuA z = (K c - K)/D = ±10~ n 
(with ± for solid and dashed lines, respectively), and n = 3, |,4, |,5 in order of decreasing T* , approaching successively 
the limit A z — [K c — K) = (dotted line). Centre: effect of including also finite left /right tunneling, A x ^ (2CK) or 
Vlr + (2IK). Shown for fixed finite 4^A Z = (K c - K)/D = ±10" 5 with 5A^/|A Z | = 2DuVlr/\K c ~ K\ = l(T m with 
m — 2, |, 1, i,0, successively approaching io-(O) = \ from above [solid lines; A 2 , (K c — K) > 0] and from below [dashed lines; 
A z , (K c — K) < 0[. Right: effect of including finite magnetic field. Shown again for 4i^A z = (K c - K)/D = ±f(T 5 , but now 
with B/4D^|A Z | = B a /\K C -K\ = f0 1/2 . As before, solid lines for A z , (K c - K ) > and dashed lines for A x , (K c - K) < 0, 
with both a =t and 4- spectra shown. Excellent agreement between NRG data and analytic curves obtained in all cases from 
a single set of fitting parameters for each model: ct = 96 and cb = 0.04 for 2CK; while ct = 0.63, cy = 2.4, cb = f .3 for 2IK. 



toward the Fermi liquid (FL) FP. Signatures of this 
crossover appear in the energy-resolved local density of 
states, since inelastic scattering ceases at energies <C T* 
where the impurity degrees of freedom are quenched. In- 
deed, the dl/dV conductance through a quantum dot 
asymmetrically coupled to source and drain leads at zero 
temperature is related Q to the scattering T-matrix: 
dl/dV oc J2a=t ii~ 7ri/lmTL,j(eV)}, where v is the lead 
density of states per spin and V the source-drain volt- 
age. The Green function is given by 

Gg£M = Gl(u)8 kk , + Gli^TMGliu), (5) 

with G;!(w) = (uj — tk + «0 + ) _1 . Conductance measure- 
ments of the 2CK device of Ref. [5| or in the proposed 2IK 
setup of Ref. [H| thus yield access to the T-matrix and 
hence the full Green function. Our goal here is to calcu- 
late the quantity t a {uS) — —TtvhnTi, a (ijj) exactly and nu- 
merically for the 2CK and 2IK models in the presence of 



symmetry-breaking perturbations described by Eqs. ([T])- 
©■ 

Results.- Preempting the technical discussion of the 
next section, we present now the key results of the Letter. 
Heralding renormalization group flow to the stable FL 
FP, the low-energy crossover scale is given generically by 

T* = c T A 2 , (6) 

where A 2 = £} =1 A? (and {A 4 ,A 5 ,A 6 } = {\%\ V B ,\* B }). 
For the 2CK model, Ai = ^A zV / 7¥, A 2 , 3 = c v vA x , y y/2K 
and A B = c B B/y/T^; while for 2IK, Ai = (K c - 
K)/VTk, A 2 + iA 3 = c v ^fL\vV LR , and X B = 
cbBs/VTk- Here ct,cv,cb = 0(1) are fitting pa- 
rameters which depend on the model and on J, and 
Tk oc is the Kondo temperatre. The perturba- 

tions associated with coupling constants A7 and Ag do 
not conserve total charge [4|, lZJ , so are ignored. 
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In the simplest case of channel anisotropy in the 2CK 
model, the result T* oc (A z ) 2 has long been estab- 
lished Q . The extension to finite Ai , A2 , A3 7^ follows by 
trivial rotation of the bare Hamiltonian in r-space, imply- 
ing directly that T* oc (A^) 2 + (A y ) 2 + (A z ) 2 , and hence 
cy = 1 for the 2CK model. However, the low-energy 
effective Hamiltonian for both 2CK and 2IK models pos- 
sesses a larger SO (8) emergent symmetry that permits 
a similar rotation, yielding the generalization, Eq. (|6|). 
Importantly, we show that the same rotation can also be 
exploited to obtain a single zero temperature Green func- 
tion for the crossover. Our exact result for the T-matrix 
is 

2muT m (u) = 1 - S aa g (w/T*) , (7) 

where G(x) = — K\ix], K[x] is the complete elliptic in- 
tegral of the first kind, yielding asymptotically Q[x] = 
1 + ix/4 - (3x/8) 2 + 0(x 3 ) for x < 1; and Q[x] = 
^(log[256x 2 ] - i^x" 1 / 2 for x > 1. S aa is the FL S- 
matrix, containing phase shift information, given by 

S 2 a c a K = (-aAi + ut\%)/\ = aSZ K , (8) 

with a — ±1 for spins t / 4- an d a = ^1 f° r channel 
L/R, such that t a (0) = h — ^ReSLa (and we have used 

Xb\\z). 

Rich physical behaviour is obtained when combina- 
tions of the relevant perturbations are applied. We now 
discuss relevant and representative cases of the general 
and exact formula, Eq. ([7]) , valid in the crossover regime 
\ui\ <C Tk- We also employ the NRG technique (for a 
review, see Ref. Q), which can be used to determine ac- 
curately t a (uj) on all energy scales, for both 2CK and 
2IK models in the presence of any perturbation. Here 
we focus on the low-energy behaviour where comparison 
can be made to the exact results. All symmetries of the 
problem are exploited, and N s = 6000 states are retained 
at each iteration. The leads of width 2D are assumed to 
have a uniform density of states v = 1/(2D), and are 
discretized logarithmically 

NRG data and exact results are presented in Fig. [1] 
with parameters given in the caption. t a (u>) vs. lo/Tk is 
plotted for the 2CK model (upper panels) and the 2IK 
model (lower panels) in the presence of various pertur- 
bations. In the left panels, finite Ai is considered (chan- 
nel asymmetry for the 2CK model, and detuning of the 
inter- impurity coupling in the 2IK model). Precisely at 
the 2CK FP, A = 0, S aa = HQ, and hence t a (w) = \ 
for |cj| -C Tk- But for Ai 7^ one immediately obtains 
\t a (u) - || ~ \u/T*\- 112 for T* < M < T K , while for 
|w| <C T*, the classic quadratic approach to the FL FP is 
given asymptotically by \t a {bS) — i<x(0)| ~ (cu/T*) 2 , with 
ta = 1 (solid lines: Ai > 0) or t a = (dashed lines: 
Ai < 0) being obtained at u = 0. This full lineshape was 
likewise obtained numerically for example in Ref. [lOj for 
the 2IK model, or for the 2CK model [IlEl] and related 



odd-impurity chains in Ref. [13|. The exact crossover 
function, Eq. ([7}, is plotted as the red line in each case, 
showing remarkable agreement for |u;| -C Tjc- 

The effect of including also left /right charge transfer 
terms (finite A2) is shown in the centre panels (Ai now 
being kept fixed). As A2 increases, the ui <C Tk line- 
shapes seen in the left panels of Fig. Q] undergo a simple 
rescaling t a (u>) — > | + |^| [t a (cj)— |] and eventually for 
I A2 1 3> I Ai| they completely flatten. In the 2CK model, 
the resulting form of t a (uj) is readily understood: rota- 
tion in r-space allows the Hamiltonian to be written in 
terms of A z only — but t a (u>) is then a weighted com- 
bination of the A z > and A z < spectra shown in 
the left panels. That the same behaviour is observed in 
the 2IK model (lower-centre panel) is a deeply non-trivial 
result however. There is no symmetry of the bare Hamil- 
tonian that permits this rotation; rather, it is the result 
of an emergent symmetry. In both cases, the crossover 
function is described perfectly by Eq. (|7J|- 

Finally, we consider the effect of applying a magnetic 
field (finite A^). In the absence of other perturbations, 
KeS ar7 = 0, and hence t^(0) = tj, (0) = \ (consistent with 
a 7I-/4 phase shift [14]). Indeed, t a (u) = | for T* < 
<C Tk since the system is near the 2CK FP. However, 
the impurity magnetization M ~ B z for the 2CK model 
(or staggered magnetization M s ~ B* in the 2IK model) ; 
thus tf((j) ^ ti(cj) since M oc /_ dui[t^(uj) — t±(u)] 7^ 0. 
For finite (staggered) magnetization therefore, the 'up' 
and 'down' spectra must deviate at finite frequency: a 
'pocket' opens between the curves at \u>\ ~ T* , whose 
area is proportional to the (staggered) magnetization. 
This behavior is observed in the right panels of Fig. [IJ 
where a fixed Ai is also included. 

For |w| -C T* , we find a linear approach to the FL 
FP, t a (u>) - t a (0) = + O , rather than the 

quadratic dependence usually associated with FL theory. 
However, this result is a perfectly natural consequence 
of the complex scattering matrix S aa . t a (uj) comprises 
contributions from both imaginary and real parts of the 
complex function Q(uj/T*) in Eq. |0, the latter of which 
contains a leading linear term. It is this contribution 
that of course dominates the low-|w| behaviour of t a (uS) 
in the presence of the magnetic field. Again, there is ex- 
cellent agreement between the NRG data and the exact 
results across the entire frequency range |w| <C Tk- We 
finally note that the conductance is obtained by averag- 
ing a =j" and \. contributions; the imaginary part of S aa 
thus cancels by Eq. ([5]). and the characteristic 'hump' 
observed in t a (ui) is absent in dl/dV. In consequence, 
dl/dV - {eV) 2 for eV < T* at finite fields. We now 
sketch the derivation of Eqs. (0 and ([SJ. 

Analytic crossover function at K 7^ K c in the 2IK 
model.- The detailed CFT analysis of the 2IK model in 
Ref. 0| demonstrated that perturbing the critical point 
by finite (K c — K) is equivalent to the action of a bound- 
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ary magnetic field h in the boundary Ising model (BIM) 
in 2 dimensions, described in the field theory limit by [f 5] 



H 



Ising 



dxe{x)id x e[x) + ihe(x = 0)a . (9) 



e(x) and a are Majorana fermion (MF) fields (a is local) 
We use the unfolded coordinate system with left moving 
convention, where x > correspond to incoming fields, 
x < to outgoing fields, and x = is the boundary itself 
(for more details see also Ref. [HI]). The renormalization 
group flow from free to fixed boundary condition [l7| is 
identical in both 2IK and BIM. The energy scale associ- 
ated with the crossover in the BIM is T* oc h 2 . We now 
apply this analogy and the machinery of CFT to obtain 
an exact result for the crossover Green function in the 
2IK model with finite (K c — K). The Fourier transform 
of the electron Green function in Eq. ([S]) factorizes into 



(l/jaaizi)^, /3 {Z 2 )) OC 



S ™' 5 «P {a( Zl )a(z 2 )) h , (10) 
(zi - z 2 )s 



where Z\ = r + ix\ and z 2 = —1x2- Here t is imaginary 
time, and x\,X2 > 0, implying that the Green function 
probes an electron propagating through the boundary. 
The electron Green function is thus given in terms of the 
two-point function of the chiral Ising spin-operator <r, 
calculated with Eq. (j9]). The magnetization in the Ising 
model evaluated at distance y from the boundary is given 
by [l8fl m(y) = {a(z 1 )a(zl)) h , where z x = iy. The full 
function m(y) at finite h has been calculated exactly by 
Chaterjee and Zamolodchikov fljl : 



m(y) oc (4 7 r/i z y)te 47r ' 1 V K {A-Kb?y) 



(11) 



where Kq is the modified Bessel function of the second 
kind. Eqs. (fTT))) , pip thus allow the Green function to 
be calculated when z 2 = z*. Extension to any z 2 7^ z* is 
possible since (<j{zi)cr(z2))h is a function of z\ — z 2 . As 
the system is not translationally invariant in space due 
to the boundary, this is a highly non-trivial result. How- 
ever, we have proved this surprising property to all orders 
in h, justifying the analytic continuation of Eq. (Ill[) to 
obtain (Vw(*i)^L(z2)) oc H m ( ~ izi + m ) ■ The 

(zi-z 2 ) 8 

special case of Eq. ([7]) with finite Ai only follows by nor- 
malization and Fourier transformation of this result. 

Generalization to arbitrary relevant perturbation. - 
Generalizing the analysis of Ref. j20j to the 2IK model 
with an arbitrary combination of perturbations {Aj}, it 
can be shown that the 2IK FP Hamiltonian becomes 
SO(8) symmetric: H = \ Y?j=i I-oo dxx 3 (x)d x Xj{x), 
where {xj} are the 8 MFs. Switching on relevant per- 
turbations at the critical point is equivalent to adding 
SHqcp = i^2j=i ^jXj(fy a > which chooses one direction 
in the 8-dimensional space. Defining a new basis in which 
only x( x ) = ^jXj{ x )/^ couples to the local opera- 
tor a, yields a Hamiltonian of the same form as Eq. (0), 



with h A and e — > x- The full Green function for 
the 2IK model with generic relevant perturbation can 
now be related to the result derived above for finite Ai 
only. The required rotation in SO (8) space is defined 
by the unitary transformation % — > Z7 xU^ = xi with 
U = cxp(j J x oo dxxi{x)x±{x)j, where 7 = arcsin^, 

A_l = \f^ 2 - Af, and x± = ^jXj/*±- Tne ke Y 

point is that the same rotation defined for the MFs can 
be used for the original fermions since linear relations ex- 
ist between their quadratic forms. The general result for 
the 2IK model is Eq. 0. 

Extension to 2CK model. - The same SO(8) represen- 
tation of the critical point is obtained in both 2IK and 
2CK models (2l|. Indeed, the CFT relevant perturba- 
tions can be matched to MFs in the 2CK model as they 
were for the 2IK model. Thus the above calculation can 
be generalized to the 2CK model, with the key results 
given in Eqs. 

Conclusion.- We have derived a single exact crossover 
Green function to describe the low-frequency crossover of 
the 2CK and 2IK models. All relevant perturbations are 
related by an emergent SO(8) symmetry, and should be 
regarded on an equal footing for T* <C Tr-, since marginal 
and irrelevant corrections to the critical point can then 
be safely neglected 22]. The derivation depends on a 



non-trivial analogy between the renormalization group 
flow in the 2IK model and in the BIM 0|, a proof of 
which is the essentially perfect agreement between the 
exact result and the full numerical solution obtained by 
NRG. 
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